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I. INTRODUCTION 

The possibility of establishing a quantum network 
of practical use is currently under active investigation. 
Some early versions of them, modest as they may be, 
have been designed and realized in real world in the re- 
cent years [H-Q , or in some instances they are under way. 
In fact, building a quantum network has been targeted 
as a fundamental goal in quantum information 0, Q an 
even a more feasible goal to accomplish than the first 
scalable quantum computer. There are other more ad- 
vanced proposals for quantum networks based on 
entanglement connections that need quantum repeaters 



11 - 13| in order to function properly and being stable 
1J-|16[. Related to this, some physical quantum mod- 
els exhibit veryremarkable long-distance entanglement 
properties |17H2(j |. Another alternative to build differ- 
ent types of quantum networks makes use of quantum 
percolation protocols [2l| - [2^ . 

Thus, it is interesting to study how different possibil- 
ities of quantum networks would behave regarding what 
we know about the world wide web. In particular, an es- 
sential ingredient in the classical network that we enjoy 
today is the ability to search web pages in the immense 
changing world that the web has come to be known. The 
key tool for performing those searches is the notion of 
the PageRank algorithm p5l - l33| . 

We envision a similar tool to perform the task of rank- 
ing quantum webpages. A quantum webpage is a node 
of the quantum network realized by means of a quan- 
tum storage device, like a quantum memory. What is 
crucially important is that a fully fledged quantum com- 
puter is not required. Thus, the quantum webpages have 
limited capabilities: storing and reading in/out quantum 
states. 

Since the notion of a quantum PageRank is by no 
means unique, it is convenient to introduce a class or 
category of possible quantum PageRanks. They must 
satisfy a set of properties that define an admissible class: 

Quantum Computer PageRank Class: 

Ql The classical PageRank must be embedded into the 
quantum class in such a way that the directed graph 
structure is preserved at the quantum level. 



Q2 The sum of all quantum PageRanks must add to 1 

Q3 The Q-PageRank admits a quantized Markov Chain 
(MC) description. 

Q4 The classical algorithm to compute the quantum 
computer PageRank belongs to the computational 
complexity class BQP. 

However, as wc have noticed that working with a real 
quantum computer seems not realistic in the near future 
and prospects for realizing quantum networks based on 
nodes made up of quantum storing devices, like quantum 
memories, are high as shown by experimental progress, 
we find also very convenient to introduce a more realistic 
class of quantum PageRanks as follows: 

Quantum PageRank Class: 

PI The classical PageRank must be embedded into the 
quantum class in such a way that the directed graph 
structure is preserved at the quantum level. 

P2 The sum of all quantum PageRanks must add to 1 

i.e. Y.^hiP^)^^■ 

P3 The Q-PageRank admits a quantized Markov Chain 
(MC) description. 

P4 The classical algorithm to compute the quantum 
PageRank belongs to the computational complex- 
ity class P. 

The crucial difference is the substitution of property 
Q4 that deals with the feasibility of calculations with 
a quantum computer in polynomial time by property P4 
which is more demanding since it means that even a clas- 
sical computer can simulate the process and find the im- 
portance of webpages in polynomial time. 

Property PI reflects the fact that originally we would 
have an internet that is a classical network represented 
by a directed graph and then shall apply some kind of 
quantization procedure in order to turn it into a quan- 
tum network. The latter must be compatible with the 
classical one, particularly preserving the directed struc- 
ture which is crucial to measure a page's authority. This 
is non-trivial and some quantization methods may fail 
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to produce a unitary quantum PagcRank importance for 
the quantum shown in Sect lllll 

With property P2 we guarantee that we have a globally 
well-defined notion of the importance of a web page at the 
quantum level. This allows us to have the probabilistic 
interpretation of the surfer's position (see Sect lllip . 

Property P3 is the key to a wide class of natural quan- 
tization methods for the classical PageRank based in the 
equivalence of this one with a classical Markov chain pro- 
cess (see Sect lIII|) . Thus, it is natural that the equivalent 
property holds true in the quantum version of the PageR- 
ank, and consequently, its description in terms of a quan- 
tized surfer's motion. 

The reason for property P4 relies on the assumption 
that we envisage a near-future scenario when a certain 
class of quantum network will be operative but not yet 
a scalable quantum computer. Therefore, we demand 
that the computation of the quantum PageRank Iq be 
efficiently carried out on a classical computer. 

In this paper we have constructed a valid quantum 
PageRank that fulfills all these requirements. We remark 
that there may be other solutions to the quantum ver- 
sion of the PageRank within the class defined above, but 
nevertheless we shall show that finding one instance of 
this quantum PageRank class is a non-trivial task. 

The definition of Class of Quantum PageRank given 
here is very general and it can accommodate very diverse 
situations, ranging from quantum channels (including en- 
tangled states) between nodes in the network to simpler 
situations where network nodes are quantum levels of a 
multilevel quantum system where information is stored 
in the same way as in the Grover algorithm. In the par- 
ticular instance of Quantum PageRank algorithm that 
we have found in subsection Results: A Quantization of 
Google PageRank , we are concerned with the latter case 
of a multilevel quantum systems where the quantum state 
is attached to a link of the network of nodes. This in- 
stance does not exhaust all possibilities represented in 
our generic Class of Quantum PageRanks by all means, 
but it is the first non-trivial instance in which sharp de- 
viations from the classical PageRank algorithm can be 
found, as explained later on in our paper. 

The explicit step-by-step description of our quantum 
PageRank algorithm is presented in SectllVI A key dis- 
tinctive feature of this quantum algorithm is that the 
importance of the quantum pages exhibit quantum fluc- 
tuations unlike its classical counterpart. These quantum 
fluctuations, as shown in the simulations in Sect|Vl show 
up in the form of time dependent importances Iq{Pi,t), 
which causes in turn that sometimes one certain pair of 
pages satisfy Iq{Pi,ti) > Iq{Pj,ti), and some other times 
the relative importance is reversed Iq{Pi, ^2) < Iq{Pji ^2), 
for time-steps such that ti < ^2- We may use an anal- 
ogy to understand this situation: the classical PageR- 
ank gives us a snapshot or photo with a fixed hierarchy 
of web pages according to their calculated importance. 
On the contrary, the quantum PageRank is more like a 
movie since the quantum importance of the pages vary 



with time. In order to produce a fixed output made of 
a list with the quantum pages sorted according to their 
importance, a natural choice we make is to compute the 
temporal average of the quantum PageRanks and their 
standard mean deviation. 

There are additional features that a certain quantum 
PageRank may have as a consequence of the definition of 
the class above. We provide hereby several useful defini- 
tions: 

Strong Hierarchical Preserving PageRank: when 
the classical hierarchical structure of a PageRank is pre- 
served upon quantization. 

This notion is too strong when the PageRank varies with 
time as it is the case at the quantum level. 
Weak Hierarchical Preserving PageRank: when 
the node with highest classical PageRank is preserved 
after quantization, but not so for the rest of pages. 
Outperforming: when the highest classical PageRank 
of a page is overcome by the quantum PageRank of that 
page. 

Outperforming may occur at one given instant or on 
average thereby leading to the natural extended concepts 
of instantaneously outperforming or average outperform- 
ing, respectively. 

In section IVII we provide a list of main results that we 
have obtained with our quantized version of the quantum 
PageRank algorithm. Remarkably, quantum fluctuations 
may change the classical hierarchy of web pages both 
instantaneously and in terms of mean values. 

This paper is organized as follows: in Sect|ll]we give an 
introduction to the classical notion of PageRank needed 
to present in Sect lllll the quantum version of it. In 
Sect llVI we present our proposal for a quantum version of 
Google PageRank algorithm. In Sect|V]wc perform nu- 
merical simulations of the quantum algorithm to repre- 
sentative directed graphs representing either an intranet 
or a general web with no special symmetry. Sect lVIl is 
devoted to conclusions. 



II. CLASSICAL PAGERANK 

Brin and Page introduced Google in 1998 [25l - [27j . a 
time when the pace at which the web was growing began 
to outstrip the ability of current search engines to yield 
useable results. A major distinction between their algo- 
rithm, called PageRank (PR) , and previous approaches is 
the fact that PR has an objective character, while other 
searchers were based up on the subjective criterium of 
the contents of the pages , because they were built as a 
collection of links that people in companies stored on a 
regular basis. In order words, PR is dynamical while the 
other approaches were static and subjective w.r.t. con- 
tents of the pages. 

The way most search engines, including Google, work 
is to continually retrieve pages from the web, index the 
words in each document, and store this information. 
Each time a user asks for a web search using a search 
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phrase, such as "search engine", the search engine de- 
termines outputs all pages on the web that contain the 
words in the searched phrase or are semantically related 
to it. 

Then, a problem arises naturally: Google now claims 
to index 50 billion pages. Roughly 95 % of the text in web 
pages is composed from a mere 10,000 words. This means 
that, for most searches, there will be a huge number of 
pages containing the words in the search phrase. What 
is needed is a mean of ranking the importance of the 
pages that fit the search criteria so that the pages can be 
sorted with the most important pages at the top of the 
list. Their success is largely due to PageRank's ability to 
rank the importance of pages in the WWW. 



A. Google PageRank 

The key idea of Google's PageRank algorithm is that 
the importance of a page is given by how many pages 
link to it. If we define I (Pi) as the importance of a page 
Pi and Bi as the set of pages linking to it, then we might 
think to put in equations the key idea put forward above 
as follows: 



HP,] 



outdcg(Pj) 



(1) 



where outdeg(Pj) is the outdegree (i.e. the number of 
outgoing links) of the page Pj. Let us define a matrix, 
called the hyperlink matrix, in which the entry in the ith 
row and jth column is: 



Hi- 



l/outdcg(Pj) 




if P, e B, 

otherwise 



(2) 



We will also form a vector / whose components are 
the PageRanks I{Pi)- The condition above defining the 
PageRank can be expressed in matrix form as: 



I = HI 



(3) 



Thus, we have recast the problem of finding the PageR- 
anks as the problem of finding the eigenvalues of a ma- 
trix [2^. We are in for a special challenge since the ma- 
trix H is a square matrix with one column for each web 
page indexed by Google. This means that H has about 
n = 50 billion columns and rows. However 7? is a sparse 
matrix, i.e. most of the entries in H are zero; in fact, 
studies show that web pages have an average of about 
10 links, meaning that, on average, all but 10 entries in 
every column are zero. 

We will choose a method known as the power method 
for finding the stationary vector I of the matrix H. How 
does the power method work? Wc begin by choosing 
a vector Jq as a candidate for / and then producing a 
sequence of vectors Ik by: 



However, as it is formulated the PageRank algorithm will 
not output a meaningful vector. We will need to patch 
the procedure in various ways. 



B. Patching the Algorithm 

It can be seen that if there are dangling nodes, pages 
that have no outlinks, then the power method will output 
the null vector. If wc consider the following example: 



whose hyperlink matrix is: 



H 





1 



and start from /q = (1, 0)* one ends up with / = (0, 0)*. 
The first patch in the tinkering of the PageRank algo- 
rithm will be replacing the column corresponding to a 
dangling node with a column of all 1 /n with n the num- 
ber of nodes. This means that virtually every dangling 
node is linking to every single node in the web, including 
itself. This prevents the power method from giving the 
null vector. This way, the disconnected graph becomes 
effectively connected at the price of giving a very low 
weight to the artificial bonds (added links). 
The graph, with the addition of the extra links would 
look like: 




with a modified hyperlink matrix, E: 



E = 



1/2 

1 1/2 



The matrix E that we obtain is, in general, (column) 
stochastic, i.e. its columns all sum up to one. From 
the theory of stochastic matrices one knows that 1 is 
always an eigenvalue. Furthermore, the convergence of 
I^ = El'^~^ to I depends on the second eigenvalue of E, 
A2. If it is smaller than 1, then the power method will 
converge. In addition, it is more rapid if IA2I is as close 
to as possible. 

Let us consider the graph in fig. [1] with E matrix: 



E = 



HP 



(4) 
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importance from Qi to G2- In order to have a meaningful 
vector J, one that has all nonzero entries, one should 
demand that the matrix be irreducible. A necessary and 
sufficient condition for it is that the graph is strongly 
connected, i.e. that given two pages there is always a 
path connecting one to the other (see [2^) chap. 8). 



Figure 1: (Color online) A graph whose matrix i?'s second 
eigenvalue, A2, is one (see text in section . 




Figure 2: (Color online) A graph that is not strongly con- 
nected, or equivalently, whose matrix E's is reducible (see text 
in section HI)) . 

One can see that the second eigenvalue of E, IA2I is equal 
to one in this case, and actually all the eigenvalues are 
on the circle of radius 1 in the complex plane. If we 
compute / with the power method starting from, say, 
Jo = (1, 0, 0, 0)* it will fail to converge. 

We will need to patch the algorithm again to ensure 
that IA2I < 1. In order to guarantee it, we will require 
E to be primitive, i.e. that there is an integer m such 
that i?™ contains all positive entries. The meaning of 
this assumption is that the graph is such that any page 
is connected by a path of at least m links to any other. 

Anticipating the interpretation of a diffusion phe- 
nomenon associated to searching the web, we can inter- 
pret the requirement of E to be primitive as the require- 
ment of finding the walker with nonzero probability on 
any site after a minimum time m. Let us now consider 
the graph in fig. [2] One can divide the graph in two 
subgraphs Gi and G2 ■ There are no links pointing from 
the subgraph G2, made of the nodes 3 and 4 to the first 
subgraph Qi, made of nodes 1 and 2. If we write down 
the matrix E: 
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1/2 


1 






we can see that there is a block that is zero, precisely the 
one that carries the information of the edges linking the 
nodes of Q2 to the nodes of Qi. 

If we calculate / with the power method starting from 
Iq = (1, 0, 0, 0)* we find: / = (0, 0, 3/5, 2/5)*. Yet, one is 
not satisfied from giving an interpretation to the vector 
/ because the nodes from the first subgraph Qi have zero 
importance albeit being linked by other nodes. This is 
caused by the reducibility of E that causes a drain of 



1. The Patched Algorithm 

In order to implement all these patches, let us imag- 
ine a walk on the graph. With probability a the walker, 
which is equivalent to the surfer mentioned earlier, will 
follow the web with stochastic matrix E and with prob- 
ability 1 — a, it will jump to any page at random. The 
matrix of this process would be: 

G:=ai?+t^l, (5) 

where 1 is a matrix with entries all set to 1 and A'^ is the 
number of nodes. This matrix is known as the Google 
Matrix. Now, the matrix G is irreducible because the 
matrix 1 is irreducible. Furthermore, it is also primitive 
since it has all positive entries. We have thus obtained a 
matrix that is both primitive and irreducible. This means 
it has a unique stationary vector that may be calculated 
with the power method. Furthermore, the result does not 
depend on the initial value Iq because the underlying 
graph is strongly connected, which is equivalent to the 
irreducibility of G, see [2^. 

The parameter a is free and needs to be tuned. It is 
known [13] that the second eigenvalue of G , A2, is such 
that IA2I < a , so one would choose a as close to zero 
as possible but in this way the structure of the web, de- 
scribed by E would not be taken into account at all. Brin 
and page chose a = 0.85 to optimize the calculation. 

2. Formulation as a Random Walk 

It is very appealing and useful to rethink the problem 
of assigning the importance of a page as the task of cal- 
culating the fraction of time a walker diffusing on the 
graph according to the stochastic process given by the 
Google matrix G. In fact we, can reformulate the Google 
PagcRank as the algorithm that computes the fraction 
of the time the walker spends on each node by defining 
the fraction on the j*'' page Tj as: 

T,=^G.,r,. (6) 

i 

Equivalently, one can say that the operational meaning of 
the PageRank algorithm is to give the probability to find 
the walker on the node Pi. Let us make it clearer by defin- 
ing a set of random variables: X^^\ X^^\ . . . , X^'^\ . . . , 
one for every time step. For each step, the random vari- 
able can take on values in the set of nodes {Pi} of the 
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web. We can recast Google PageRank in the language of 
a Markov Chain. Thus, from eq. ^ written as: 

Pr(X("+i) =P.)=J2 G«,Pr(^'"' = Pj) (7) 
j 

and from the law of total probability: 

Pr(X("+i) = PO = 
2^.Pr(X("+i) = Pi|X(") = P,)Pr(X(") = Pj), 

one can interpret the stochastic matrix G as the condi- 
tional probability linking one time step to the other, i.e.: 

Gij = Pr(X("+i' = P,|X("' = Pj). (9) 

We will make use in the following of the latter interpreta- 
tion of Google PageRank to devise methods to quantize 
it. 



III. QUANTUM PAGERANKS 

Quantum walks in their discrete time formulation were 
known already to FeynmanJSJ] and since then, they were 
rediscovered many times [35| and in contexts as different 
as quantum cellular automata [IB, US] and the halting 
problem of the quantum Turing machine [H, [s^ . For 
simplicity, let us discuss possible ways for quantizing a 
quantum walk taking place on the line. Later, we shall 
generalize it to an arbitrary graph. From now on, we 
shall make the notation lighter denoting each node (or 
page) Pi simply by i as shown in the following figure: 




The naive way of quantizing this random walk would be 
to go from the index set {i\i G Z} to a Hilbert space of 
states H — span{| i ) | i G Z} as shown in figure [H 
Following the key idea outlined above, one could define 
the quantum importance of a page as: 

/(P)=Pr(X = P):=|(z|V)r (10) 




Figure 3: (Color online) The naive way to quantize a random 
walk on a line. 



where | ■0 ) is the state of the system after it has diffused 
on the graph. However, this quantization procedure is 
not viable. This is because the direct quantization of the 
time step evolution operator as 

u ^^\i + i){i\ + ^/T~l^\i~i){i\ (11) 

is not unitary. Indeed, while (0|2) — one has that 
(0|[/tC/|2) 7^ in the general case of p e [0, 1]. 
This difficulty can be overcome enlarging the Hilbert 
space. There are various ways to do it: 

• adding a coin space "He to each site on the quantum 
network. A coin space encodes the possibility to go 
left or right i.e. He = span{|L), |R)}. 

• defining the Hilbert space as the space of (ori- 
ented) edges and treating the vertices as scattering 
centers. The resulting walk is called a Scattering 
Quantum Walk. 

• using Szegedy's [ioj procedure to quantize Markov 
chains. 

We will discuss the latter way for it will give us a 
valid quantization of Google's PageRank that satisfies 
the properties 1 — 4 introduced in SectU 

3. Szegedy 's Quantization of Markov Chains 

We have seen that the Google matrix G can be seen 
as the time step evolution operator of a Markov chain, 
or equivalently, of a discrete-time classical random walk: 
the two terms will be used interchangeably from now 
on. Szegedy put forward a general scheme to quantize a 
Markov chain. Let G be a x iV stochastic matrix rep- 
resenting a Markov chain on an A- vertex graph. In order 
to introduce a discrete-time quantum walk on the same 
graph we use as the Hilbert space the span of all vectors 
representing the N x N (directed) edges of the graphs i.e. 
n = span{\i) i\ j) 2, with i,j e Af x A} = ® C^. The 
order of the spaces in the tensor product is important 
here because we are dealing with a directed graph. We 
stress this fact in the notation using subindices 1 and 2. 
Let us define the vectors: 

N 
k=l 

that is a superposition of the vectors representing the 
edges outgoing from the j'*'' vertex. The weights are given 
by the (square root of the) stochastic matrix G. 

One can easily verify that due to the stochasticity of 
G the vectors for j = 1, 2, . . . , N are normalized. The 
operator 

JV 

n:=E|V',)(V,|, (13) 
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is then a projector onto the subspace generated by the 
vectors j^/'j) for j = 1, 2, . . . , N. With this, the step of the 
quantum walk is then given by 

U := S{2U-1) (14) 

where S is the swap operator i.e. 

N 

S=Y.\j,k){k,j\. (15) 

Before continuing let us point out that the time step 
operator is unitary: 

uu^ = S'(2n-i)(2n-i)S''f = 5(4n2-2n-2n+i)S'^ = i 

from the unitarity of S and the fact that 11 is a projector 
and squares to itself. Analogously: 

u^u = (2n - i)S'^s'(2n - 1) = 4^ - 2n - 2n + 1 = i 

(17) 

The time step is thus the effect of a coin flip followed by 
a swap operator. Let us look more closely at the coin flip 
operation: 

" / 1 \ 

2n-l = ^ 2|V',)(^',|-^l). (18) 

The vectors contain the information of the di- 
rected links that connect the j*'' node to all its neighbors 
to which it is connected through the stochastic matrix G. 
The sum on all nodes of the operators 2 \ipj){ipj \ — is 
nothing but a reflection around the subspace spanned by 
the vectors jV'j) and has the effect of enhancing the am- 
plitudes of the mentioned directed edges at the expenses 
of the others. Furthermore, the swap operator preserves 
the unitarity of the step. 

4- Solving the Eigenvalue Problem for the Walk Operator 

The quantization procedure based on the unitary op- 
erator (fT4)) allows us to get remarkable insight on the 
properties of the walk by a systematic analysis of the 
spectral properties of U. The spectrum of the quantized 
walk is related to the spectrum of the original stochastic 
matrix that in turn we have seen has a key role on the 
properties of the related classical random walk. 

Let us define the following N x N matrix D that will 
play an important role in relating the spectra of the clas- 
sical and quantum walks by specifying its entries as fol- 
lows: 

Dij := ^y Gij Gji , (19) 

where there is no sum over repeated indexes. Let us also 
define the following operator A from the space of vertices, 
C^, to the space of edges, C^: 

N 

A:=Y.\i:,){j\ (20) 



It has the following properties, which are straightforward 
to prove: 

1. A'fA = 1 

2. AA^ = n 

3. A^SA = D 

The matrix D is symmetric by construction. The eigen- 
value problem D\X) = A|A) can in principle be solved 
yielding N eigenvalues A with the associated eigenvec- 
tors I A): 

a{D):^{X,\X)}. (21) 

Let us now consider the the following vectors out of 
them: |A) := A\X) on the space of the quantized Markov 
chain i.e. C'^. In order to obtain the spectrum of 
U, we will first isolate an invariant subspace for U and 
look for eigenvalues and eigenvectors in this space. We 
will then concentrate on the orthogonal complement of 
it. We will argue that the interesting part of the Hilbert 
space for the dynamics is the aforementioned invariant 
subspace. In order to identify the invariant subspace of 
the walk operator let us see the effect of U on |A): 

C/|A) = S'(2AAf - 1)A|A) = 5|A), (22) 

where we used the property 1 and 2. Let us see also its 
effect on S\X): 

US\X) = S{2AA^ - 1)SA\X) = (2AS' - 1)|A), (23) 

using properties 2 and 3 and the fact that the vectors |A) 
are eigenvectors of D. From (1^^ and we can deduce 
that the subspace 

Iu:^{\X),S\X)}, (24) 

is invariant under the walk operator U. It is thus sensible 
to solve the eigenvalue problem: 

u\fl)^^i\^l), (25) 

for the walk operator restricted to the invariant subspace 
Following what we have said, let us make an edu- 
cated ansatz for the eigenstates of U: 

|m) = |A)-m^|A). (26) 

We have that: 

u\^i)^^i\x) + {l-2^lX)s\x), (27) 

thereby the condition for |/.() to be eigenstate of U is 

^fi' = {l- 2fiX), (28) 

which yields /i = A ± VA^ — 1 = exp(±i arccos A). 

We note also that the span of the vectors |A) coin- 
cides with the span of the vectors \tpj). Indeed we have 

Ea \~^)CM = ^Ea |A)(A|At = n = mi^jl ■ 
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To complete our analysis let us point out that on the 
orthogonal complement to the span of the vectors {ipj), 
the action of the walk operator U = 5(211 — 1) is just 
—S, which has eigenvalues ±1. This is because 11 yield 
the null vector when applied to vectors belonging to this 
subspace. We conclude that the spectrum of U is the set 



a{U) {±1, exp(±i arccos A)}, 



(29) 



where A are the eigenvalues of D. In the following we 
will need the quantum walk where two steps at a time 
arc performed, with operator U^. We can advance that 
the interesting subspace. where we have dynamics, is the 
span of the vectors and Sltpj) where the walk oper- 
ator acts nontrivially. 

Furthermore, considering the two-step evolution opera- 
tor, one can see that: 



= (2S'ns'- i)(2n- 1). 



(30) 



Therefore under the two-steps operator the Hilbert 
space splits naturally into the subspaces Hdyn = 
span{|-0j), Slipj)} where dynamics takes place and its or- 
thogonal complement T-Lnodyn = T^dyn '"^^ seen 
from ((30|) C/^ acts trivially in such a way that, obvi- 
ously, T-L = T-Ldyn © T-Lnodyn- The dimension of T-Ldyn 
is at most 2N and, remembering that the dimension of 
T-i = (g) is -/V^ we can conclude that the spectrum 
of corresponding to Hdyn is composed by, at most, 
the 2N values 



{exp(±2i arccos A)}, 



(31) 



and the rest of the spectrum, corresponding to T-Lnodyn 
where acts trivially, is composed of at least A^^ — 2N 
I's. 

The analysis presented above will allow us to save 
computational resources when calculating the Quantum 
PageRank because of the presence of the invariant sub- 
space T-Ldyn- A different type of problem is concerned 
with using a quantum computer to perform a quan- 
tum computation that might improve the efficiency of 
the classical PageRank algorithm. An adiabatic quan- 
tum computation can be set to to compute the classical 
PageRank vector. In this case the classical PageRank 
vector is encoded in a quantum state and an adiabatic 
quantum computation is described in order to reach this 
state m. 



IV. A QUANTIZATION OF GOOGLE 
PAGERANK 



In this section we define a valid Quantum PageRank 
and take advantage of the analysis presented above to 
provide an efficient algorithm for its calculation as re- 
quested in Sectm A natural way to define a quantiza- 
tion of the importance of a node or page in the quantum 
network associated to a directed graph is to exploit the 



connection with the Markov chain process in which the 
fictitious walker is now subjected to quantum superpo- 
sitions of paths throughout the quantum web. In this 
way, the instantaneous importance of a quantum web 
page, denoted as Iq{Pi, m), is given by the probability of 
finding the walker in that page Pi at the node i of the 
network after m time steps. As we have said before, the 
Hilbert space of this quantum walk is the set of directed 
links of the graph, Ti = (E) where the numbering of 
the vector spaces is meaningful due to the directedness of 
the underlying graph and the second space in the tensor 
product contains the information of where the directed 
link points to. It seems natural then to project onto a 
vector of this second space \i)2 obtaining the quantum 
state \Iq{Pi,m)) that contains the superposition of the 
nodes that were linked to it. To quantify the importance 
one can then extract a positive number calculating the 
norm of \Iq{Pi, to)), and with this we obtain an instanta- 
neous list of page ranks including quantum fluctuations 
of the network. Thus we expect its instantaneous value 
to oscillate in time as a result of the underlying coherent 
dynamics. 

The method for computing the instantaneous PageR- 
ank of the page Pi is to start from an initial vector |'(/'(0)) 
and to let it evolve according to the two-step evolution 
operator (in order to swap the directions of the edges 
an even number of times, thus preserving the graph's 
directedness). Then, we need to project onto \i)2, and 
finally to take the norm of the resulting quantum state: 

Iq{P„m) = {m\ C/^""N)2(*|C/'"|V'(0)). (32) 

In order to implement the full procedure, one starts from 
the stochastic matrix G representing the classical Google 
search that we want to quantize, forms the matrix D and 
obtains its spectrum (t{D) = {A}. One then forms the 
states I A) = ^|A), in terms of which the eigenvectors of 
the walk operator, in the subspace where the dynam- 
ics takes place are written. 

Using the spectral decomposition of U one can then 
arrive at a closed analytical expression for our quantum 
instantaneous PageRank: 



(33) 



where = |A) — /i^jA), fj, = exp(±i arccos A) and 
I ■0(0)) is taken to he on the dynamical subspace. Note 
that due to the fact that \i)2 are a basis of C^, in other 
words J2i N)2(«| = 1 , from ((5^ one can see: 



(^(o)|t/t'™IEJ*)2(*|c/'™IV'(o)) = 1 Vto, 



(34) 



meaning that in the quantum version of Google PageR- 
ank we have that the normalization condition (|34p is 
preserved at all times allowing to interpret the quantity 
Iq{Pi, to) as the instantaneous relative importance of the 
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page Pi, thereby reproducing a basic sum rule that also 
holds in the classical domain. 

In order to integrate out the fluctuations arising from 
the coherent evolution we also introduce the average im- 
portance of the page i sitting on the i*'' node {Iq{Pi)) 
as: 

■■=mY. w^^)- (35) 

m=0 

We also use its variance or standard mean deviation: 

AI,{P,) := y^(/2(i^0> - iW))', (36) 
as a measure of its fluctuations. 

In order to obtain the Quantum PageRank values of the 
nodes of a digraph we apply the algorithm that comes 
out of the analysis presented above. Namely the steps 
one has to perform are: 
Quantum PageRank Protocol 

Step 1/ Write the Google matrix for the digraph Q. 

Step 2/ Write down the matrix D (see eq. (|19p) and 
calculate its eigenvalues and eigenvectors. 

Step 3/ Find the eigenvectors and eigenvalues of the 
two-step quantum diffusion operator in the dy- 
namical subspace Hdyn (see section [ni] for the de- 
tails). 

Step 4/ Extract the Quantum PageRank value in 
time (|33|1 . its mean ([35|) and standard deviation 
P6p starting from the initial condition — 



V. RESULTS: SIMULATIONS FOR QUANTUM 
PAGERANKS 

After developing a quantum version of the Google 
PageRank in the previous section, it is necessary to apply 
it to specific networks and by means of simulations, see 
how it behaves as compared with the classical algorithm 
of PageRank. 

We have put to test our new quantum version of the 
PageRank algorithm in the case of a binary directed tree 
with 3 levels (see fig. |4]) and a small albeit general, with 
no special property, directed graph (see fig. E]). We will 
describe the results in subsections IV Al and IV Bl respec- 
tively. 

A. Case Study 1: A Tree Graph 

In this subsection we will display the results obtained 
in the case of a tree graph (see fig. H]). This type of 
directed graph has a clear meaning in terms of a web 
network: it represents an intranet with the root node 




Figure 4: (Color online)The three level tree considered in 
the text to benchmark the Quantum PageRank (see text in 
section |V| . Each node represents a web page in an intranet 
with the root node being its home page. 




Figure 5: (Color online) The general graph with 7 nodes con- 
sidered in the text for benchmarking the Quantum PageRank 
(see text in section |V]) . Each node represents a web page in 
this directed quantum network. 



being the home page of a certain website and its leaves 
representing internal web pages. This case of study has 
been extensively studied classically [Slj. The quantum 
algorithm presented above is implemented numerically. 
The quantum PageRank of the root page clearly oscil- 
lates in time and attains values that are higher than the 
classical counterpart, as shown in Fig. [6l According to 
the properties studied in Sect HI our quantum PageRank 
is instantaneously outperforming. 

It is rather remarkable that a quantum version of the 
PageRank may have an enhancement of the importance 
in the home page with respect to its classical counter- 
part. This is achieved merely by quantum means, with- 
out changing the connectivity of the original directed 
graph as has been proposed classically [3l|. 

As for the quantum fluctuations present in the im- 
portance of the root page, it is important to emphasize 
that they remain bounded during the evolution as can be 
checked from Fig. [SI In addition, the classical value is al- 
ways inside the range of the quantum fluctuations. This 
feature is found to be true for all pages in the directed bi- 
nary tree network conflguration (see FigjT]). A distinctive 
feature of the Quantum PageRank is that the hierarchy 
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Figure 6; (Color online) The evolution of the instantaneous 
Quantum PageRank Iq with time for the root page (home 
page) in the case of a directed binary tree in Fig|4l 
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Time 

Figure 7: (Color online) The evolution with time of the instan- 
taneous Quantum PageRank algorithm Iq defined in Sect llVI 
for web pages (nodes) in the case of a directed binary tree 
graph in FigU Only one page per level of the tree is dis- 
played because pages that are in the same level have equal 
Quantum PageRank. 

is not preserved at every time. Figure [7| clearly display 
the crossings of the instantaneous Quantum PageRanks. 
In order to extract a fixed value for the relative impor- 
tances of the pages, we compute the mean value of the 
instantaneous quantum PageRank in time from eq. (|35|) 
and its variance from eq. p6p . 

One can notice that the hierarchy of the pages is pre- 
served on average and that the errors i.e. the variances 
are negligible compared to the means. We can thus infer 
the pages' hierarchy as predicted by the classical algo- 
rithm from the Quantum PageRanks' averages as can be 
clearly seen in figure [8] for this particular case of directed 
binary tree graph. 

Furthermore we can calculate a coarse grained evolu- 
tion in time of the instantaneous Quantum PageRank. 
We divide its total evolution time of M steps in L equal 
segments made of an integer number M/L steps. We 




0.05 ' ' ' ' ' ' ' ' ' 

012345673 

Nodes 



Figure 8: (Color online) Comparison of the hierarchies that 
result from the Classical and the Quantum PageRank in the 
case of the directed graph tree shown in FigU) Error bars 
in the quantum case are computed with the standard mean 
deviation ((35)l . (l36ll . 



Level 


Classical PageRank 


Average Q- PageRank 


Variance 


1 


0.37291 


0.355905 


0.0156461 


2 


0.18012 


0.151437 


0.0067747 


3 


0.06671 


0.085305 


0.0022797 



Table I: Comparison Classical and Averaged Q-PageRank in 
the case of the tree graph shown in Fig. [l] Given the symme- 
try of the tree only the values for each level are displayed. 



calculate the mean in every segment: 

^ (n+l)M/L-l 

J,(P„n):=— W,m) ,n^O,...,L-l. 

ni—nAI/L 

(37) 

The result of integrating out the oscillations in a coarse 
grained time step is shown in figure [9l An interesting 
feature of this case is that the quantum PageRank oscil- 
lations shows a modulation with a clearly visible envelope 
for each node. It strongly enforces the idea that the pro- 
posed algorithm has a solid grounding on the classical 
PageRank algorithm and it is a valid quantization of it. 



B. Case Study 2: A General Graph 

We have performed a calculation of the quantum 
PageRank also in the case of a general directed graph 
with no particular symmetry (see fig. [5|). 
The value of the instantaneous quantum PageRanks are 
found to display oscillations that are bounded and in- 
clude the value found by the classical PageRank. The 
QPR of the node with the highest classical PageRank at- 
tains, at given times, values of the QPR that are higher 
than the classical counterpart (see fig. [TO]). As seen in 
the case of the binary tree treated above the Quantum 
PageRank is found to be instantaneously outperforming 
according to the definition given in section |T1 
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Coarse grained time 

Figure 9; (Color online) The evolution in a coarse grained 
time (see text in section [V]! of the Quantum PageRanks in 
the case of the directed graph tree shown in FigU Again 
only one page per level is being displayed because pages that 
are in the same level have equal Quantum PageRank. 
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Time 



Figure 10: (Color online)The evolution in time of the Quan- 
tum PageRank Iq of pages 7, 5 (that classically have the high- 
est and nearly degenerate PageRank) and of page 4 (that clas- 
sically has the lowest PageRank) in the case of general graph 
shown in FigO 



The classical hierarchy is not preserved by the QPR at 
any given time (as is clearly shown by Fig. I10[ see cap- 
tion) . We can notice crossings in the importance given by 
the instantaneous Quantum PageRank even between the 
pages that classically have highest and lowest PageRank. 
Furthermore, the nodes that have a very close classical 
PageRank are shown to have a very similar behaviour in 
time of their instantaneous Quantum PageRank. 

Remarkably enough, we find that the Quantum PageR- 
anks' averages do not give us the same hierarchy as in the 
classical case (see Fig. [TT|) . Nevertheless, it is possible 
to clearly distinguish, within the error bar given by the 
variance, which pages have highest and lowest classical 
PageRank. 

The analysis with a coarse graining in time of the in- 
stantaneous Quantum PageRank (see Fig. [12]) reinforces 
the conclusion that the classical PageRanks are still dis- 
tillablc in the case of pages with highest and lowest clas- 
sical importances. 



Node 


Classical PageRank 


Average Q-PageRank 


Variance 


1 


0.051019 


0.089076 


0.0021759 


2 


0.061860 


0.126546 


0.0050376 


3 


0.077924 


0.130587 


0.0040337 


4 


0.028940 


0.076586 


0.0014675 


5 


0.362387 


0.217691 


0.0111097 


6 


0.047981 


0.131345 


0.0049477 


7 


0.369889 


0.228169 


0.010549 



Table II: Comparison Classical and Averaged Q-PageRank in 
the case of the general graph shown in fig. [S] 




Nodes 



Figure 11: (Color online) Comparison of the hierarchies that 
result from the Classical and the Quantum PageRank in the 
case of the general graph shown in Fig[5l Error bars in the 
quantum case are computed with the standard mean devia- 
tion (I35l).(l36l). 
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Figure 12: (Color online) The evolution in a coarse grained 
time (see text in section |V]) of the Quantum PageRanks in 
the case of the general graph shown in FigO 
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VI. CONCLUSIONS 

In this paper we have proposed a notion of a class 
of protocols that qualify to be considered a quantum 
version of the classical PageRank algorithm employed 
by the Google search engine (see Scct|T|. In addition, 
we have constructed a step-by-step protocol explicitly in 
Sect llVI based on the quantization of Markov chains for 
directed graphs. This is a non-trivial problem since deal- 
ing with quantum versions of digraphs may produce uni- 
tarity problems (see Scct lIII|) . We have tested our quan- 
tum PageRank algorithm with two web networks in order 
to gain insight into the specific behaviour of this protocol. 
One network is a binary tree graph representing an in- 
tranet with the root being the home page of it. The other 
is a general directed graph with no specific structure. 

From our numerical simulations we have found that 
our quantum PageRank has very interesting properties. 
For the directed binary tree: 

i/ The quantum PageRank for the root page is instanta- 
neously outperforming with respect to the classical value. 
This is a manifestation of the quantum fluctuations in- 
herent to the quantum version of the algorithm and al- 
lows us to have an enhancement of the importance of the 
root page without changing the topology of the original 
network. 

ii/ The instantaneous values of the quantum PageRanks 
for the nodes violate the hierarchy of the classical values, 
iii/ The mean values of the quantum PageRanks includ- 
ing its standard deviation preserve the hierarchy of the 
classical values. 

For the general directed graph: 

i/ The quantum PageRank for the web page with the 
highest classical PageRank is some times higher than the 
classical values obtained with the standard algorithm. 
This means that the quantum version of the PageRank 
is instantaneously outperforming with respect to the clas- 
sical value. 



ii/ The instantaneous values of the quantum PageRanks 
for the nodes violate the hierarchy of the classical values, 
iii/ Remarkably enough, there are pages with mean val- 
ues, including its standard deviation, of their quantum 
PageRank that violate the hierarchy of the classical val- 
ues. 

These properties are a clear manifestation that our pro- 
posal for a quantum version of the PageRank algorithm 
exhibits nontrivial features with respect to the Classical 
PageRank. 

As the main purpose of our work is to devise a quantum 
PageRank algorithm by overcoming certain difficulties 
explained in the paper, thus far we have dealt only with 
small networks representing different types of web. It 
would be very interesting to perform computations with 
the quantum PageRank applied to very large networks 
with the properties exhibited by the complex structure 
of the real web 

An interesting issue is whether the classically first ranked 
page remains with the highest importance also at the 
quantum level. While we have shown that generically 
this is not the case for instantaneous values of the quan- 
tum PageRank (33) in Fig. 10 and Fig. 13, they remain 
first-ranked with the time-average values of the quantum 
PageRank (35). It remains open to see what happens 
for larger networks having very similar first-ranked nodes 
when the effect of quantum fluctuations are taken into ac- 
count. This will depend on the topology of the lattice as 
well. 
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